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It is commonly believed that the observed solar activity is driven by the dissipation of free 
(nonpotential) magnetic energy injected into the corona by dynamic processes in the photosphere. 
The enormous range of scales involved in the solar activity makes it difficult to track down 
the photospheric origin of each coronal dissipation event, especially in the presence of complex 
magnetic topologies. In this paper, we propose a new statistical-physical approach for testing 
the photosphere-corona coupling as manifested in a quiet solar region. We investigate large 
ensembles of photospheric and coronal events detected in co-aligned sets of images provided 
respectively by SOHO MDI and STEREO EUVI instruments. We show that for properly adjusted 
detection thresholds corresponding to the same degree of intermittency in the photosphere and 
corona, the detected events are described by similar occurrence probability distributions of time- 
integrated quantities but significantly different geometric properties. We derive a set of scaling 
relations reconciling these empirical results and enabling statistical forecast of coronal dynamics 
based on photospheric measurements. The perform analysis suggests that multiscale energy 
dissipation in the corona is directly controlled by the turbulent photospheric convection so that 
probabilistic properties of energy release events in the photosphere are imprinted in the corona. 
The nonlocal nature of the photospheric network makes this coupling essentially nonlocal and 
non-deterministic. Our results are in agreement with the Parker's coupling scenario in which 
random photospheric shuffling generates marginally stable magnetic discontinuities at the coronal 
level, and are also consistent with an impulsive wave heating involving multiscale Alfvenic wave 
packets. More research is needed to tell the two mechanisms apart. 
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ABSTRACT 



The Sun is an inherently complex and mul- 
tiscale system which continues to challenge the- 



with scales typically on the order of a solar radius 
are the sources of high speed solar wind streams. 
Coronal active region magnetic loops, with char- 
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acteristic scales of 10,000 km, exhibit somewhat 
steady heating. Direct current heating of the coro- 
nal plasma, unresolved with present instruments, 
takes place at a typical ohmic dissipation scale in 
the corona, which is likely to be as small as 100 km. 
The Sun also exhibits a broad range of temporal 
scales. The largest features are observed to last 
for months. Coronal holes and coronal streamers 
can persist for several 27-day solar rotations with 
little change, while M- and X-class solar flares re- 
lease a vast amount of energy essentially doubling 
the solar lumino sity for a brief period of a few 
minu tes (see e.g. Aschwandenll2006 : Solanki et al 
20061 and references therein). 



Despite the enormous complexity of solar pho- 
tosphere and corona, the fundamental underly- 
ing process controlling their interaction could be 
surprisingly simple. It is widely believed that 
the ultimate power source for all energy release 
processes in the upper solar atmosphere is the 
convection at the solar surface mediated by mag- 
netic field. The observed solar activity is likely 
to be driven by the dissipation of free (nonpo- 
tential) magnetic energy which is injected into 
the corona d ue to a variety of processes in the 
photosphere (jKlimchukl 120061 ) . Among the most 
important such processes are the injection of hc- 
licity associated with new magnetic structures 
dFalconer et al.l 120081: lAbramenko fc Yurchvshvnl 
l2O10HMcAteer et al.l2010UConlon et al.ll2010l ) and 
the fragmentation of the existing magnetic flux 
subject to magnetic footpoint shufflin g due to the 



turbu l ent photospheric con v ection (ILa mb et al 
l2008t IChaouche etaH 1201 it lAbramenkol 120051 ) 



The amount of free coronal magnetic energy sup- 
plied through these processes seems to be in a 
quasi-steady balance with the energy released via 



coron al heating and flaring activity (ISimon et al 



200 If ). 



Although the convective scenario of the pho- 
tosphere - corona coupling has been reproduced 
in numerical simulations, its empirical verifica- 
tion presents a significant challenge. The broad 
ranges of scales involved in solar activity makes 
it impossible to track the photospheric origin 
of each particular coronal dissipation event, es- 
pecially of the events associated wit h complex 
multisclae magnetic topologies (e.g., Antiochosl 



19981: lAntiochos et al.lf2007tlschriiver et all 1 199a 



2002[lGeorgoulisl2005HMcAteer et "all2010t|Parnell et al 



20091 : lUritskv fc Davilal l2012h . Coronal bright- 
enings occurring on strongly entangled mag- 
netic field lines tend to interact in a nonlinear 
fashion producing vast regions of secondary in- 
stabi lities in respon s e to relatively small trig - 
gers (|Klimchukll2f)0i Fucntcs k Klimchuljl2010h . 
The bulk of the emitting coronal plasma is not 
steady b ut is rather dynamic and co nstantly 



evolving (|Viall fc Klimchukl l201li 120121) . Tur 



bulcnt behavior of high- Reynolds num ber pho 



tospheric flows (jAbramenko et al.l 120021 ) as well 



Balke et al.lll993 : Lawrence et al.ll993t Abramenko et al 



as turbulence and fract al diffusion in the corona 
dBerghmans et al.lll998t lUritskv et "aDl2007l [2009t 
lAschwandenl l2012al)~ further increases the complex- 
ity of the coupling between the two solar regions. 
In addition, a noticeable fraction of magnetic dis- 
continuities created by the photospheric motion is 
likely to remain silent until a local plasma insta- 
bility condition is met ( Vlahos fc Georgoulisll2004 ; 
lAschwandenll2006[ ). 

Based on this picture, the coronal response to 
a single reconfiguration of the photospheric field 
is likely to be nonlocal, non-simultaneous, and en- 
ergetically disproportional to the triggering pho- 
tospheric event when investigated on a case-by- 
case basis. Searching for driving signatures of in- 
dividual coronal brightcnings in relatively station- 
ary quiet solar regions can be particularly difficult, 
while finding deterministic cause and effect links 
governing an ensemble of such events is close to 
impossible. 

In this paper, we adopt a different strategy fo- 
cusing on statistical-physical signatures of photo- 
spheric an coronal events occurring in a quiet Sun. 
We consider two types of events - photospheric 
magnetic field enhancements and coronal heating 
events observed simultaneously within the same 
solar region. The events are detected by threshold- 
ing SOHO MDI line of sight photospheric magne- 
tograms and STEREO EUVI coronal emission im- 
ages, correspondingly. The events represent inter- 
mittent conditions in the photosphere and corona, 
with the net area covered by the events at any 
given time being less then five percent of the total 
image area. 

The results of our analysis show that for the de- 
tection thresholds leading to same degree of inter- 
mittency of the photospheric and coronal events, 
they are described by nearly identical power-law 
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probability distributions of event areas, spatiotem- 
poral volumes, and integrated intensities. The 
power-law distribution exponents characterizing 
solar photosphere and corona undergo similar vari- 
ations during the studied 14-hour interval. The 
statistical signatures obtained suggest that tem- 
poral and spatial scales of coronal heating events 
can be directly controlled by convective shuffling 
motions of the photospheric magnetic field. The 
results are in agreement with the Parker's cou- 
pling scenario in which multiscale photospheric 
flows generate a large amount of marginally stable 
magnetic structures at the coronal level, the free 
magnetic energy stored by the structures being to 
a lar ge extent predefined by the photospheric mo- 



tion (jParkerlll983l Il988l ). The coronal brighten- 



ings release the energy of a random subset of these 
magnetic discontinuities repeating their proba- 
bilistic pattern. Alternatively, the results can in- 
dicate that the photospheric dynamics generates 
a broad-band spectrum of Alfven waves trans- 



mitted into the chromos 


ahere and eventually into 


the corona (e.g., Davilal 


1987: Ofman et al.l 1998: 


Kudoh & Shibata! 


119991: 


iDe Pontieu et al.l l2007h. 



In either case, our observations suggest strong en- 
ergy coupling between ensemble-based properties 
of solar photosphere and corona - the effect we 
refer to as stochastic coupling. 

2. Methods and results 

We have studied simultaneously collected sets 
of SOHO MDI line o f sight (LOS) magnetograms 
dScherrer et al.lll995l) and STEREO-B EUVI im- 
ages ( Howard et al. 20081 ) showing the dynam- 
ics of a quiet Sun during 17:29:00 04/05/2007 - 
10:58 04/06/2007. The images were preprocessed 
and dcrotatcd using the standard Solar Soft soft- 
ware. STEREO EUVI images were rebinned to 
down to the spatial resolution of SOHO MDI 
data (0.6 arcsec, or about 0.44 Mm) and inter- 
polated in time using linear interpolation scheme 
to match the timing of the nearest MDI snapshot. 
To reduce the influence of the granulat ion noise 
(|DeForest et al.l 120071 : lLamb et all 120081 ) the im- 
ages were also smoothed in time using a 5 minute 
boxcar moving window. The obtained data cubes 
contained 770 x 500 x 952 data points in the lati- 
tudinal, longitudinal and temporal directions, cor- 
respondingly. As a result of the described co- 
alignment procedures, both data sets (MDI and 



EUVI) were described by the same sets of helio- 
graphic positions and times. 

Fig. Q]shows full-disk SOHO MDI and STEREO 
EUVI images with the studied field of view (FoV) 
outlined by the white frame. The overall solar 
activity level remained very low, with the GOES 
X-ray fluxes staying below 10~ 8 W m~2. Two 
decaying active regions, NOAA 10949 and NOAA 
10950, were seen outside of the FoV; neither of 
them produced solar flares starting from several 
days before till several days after the studied time 
interval. The position of the FoV close to the disk 
center ensures small projection distortions. The 
latitudinal size of the region (about one fourth of 
the visible solar diameter) is sufficiently large to 
obtained statistically representative measures of 
the photosphere - corona interaction. 

The magnetic network and coronal emission 
pattern of the studied region are typical of a 
quiet Sun (Fig. [2]). The LOS magnetic field 
features a hierarchy of spatial scales representing 
several levels of photospheric granulation with a 
roughly uniform distribution of large-scale mag- 
netic flux across the FoV. The imbalance between 
positive and negative fluxes is about 15% of the 
total unsigned flux in this region, with the preva- 
lence of the negative flux. The EUVI emission 
field exhibits a diffuse, fractal-like spatial distri- 
bution of intensity characteristic of a nonflaring 
corona. Some of the bright coronal regions co- 
incide with mixed-polarity magnetic regions, in 
agreement with earlier stu dies of the quiet Su n 
magnetic network (see e.g. iFalconer et al. 1998 ). 
The appearance of bright coronal points in the 
vicinity of magnetic field reversals can be a sig- 
nature of low-altitude magnetic reconnection oc- 
curring on magnetic disc ontinuities driven by the 
photospheric motion (see Aschwandenll 20061 chap. 
10 for a review). 

To study the dynamics of photospheric mag- 
netic field and coronal emission flux, we applied 
the same numerical approach which enabled a 
detection of spatiotemporal intermittent events 
in each data sets. The approac h is based on 



spatiotempo ral tracking algorithm (jUritskv et al 



20071 l2010al) which treats photospheric or coronal 
features staying for more than one image sampling 
interval above a specified detection threshold as 
subvolumes in the three-dimensional space-time. 
The time history of each feature is represented by 
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spatially and temporally connected sets of image 
pixels. 

The idea of the threshold-based spatiotemporal 
feature tracking has first been imp lemented in the 
contex t of solar magnetogr ams by iBerger fc Title 



(|l996h : IBerger et all (jl998l ). The object tracking 
method used in their study was based on mea- 
suring the centroid locations of the most compact 
and intense magnetic features. Here we invoke 
a different approach which does not rely on the 
centroid location as a guiding center for feature 
tracking. Time-evolving magnetic concentrations 
are treated as true 3-dimensional objects, which 
makes it possible to study a more general class 
of magnetic and emission structures whose posi- 
tions cannot be adequately described by a cen- 
troid. Compared to more recent works, our al- 
gorithmic solution reverses the usual order of fea- 
ture tracking steps (feature segmentation within in 
an individual image followed by cross-frame fea- 
ture identificat i on wi thin an image set, see e.g. 
DeForest et al. ( 2007t )V Instead, we start with 
detecting temporal traces of features in individ- 
ual solar locations, and then consider their spatial 
adjacency. We found this algorithmic solution to 
be remarkably fast and memory-efficient, with the 
size of the output data array being typically at 
least by two order of magnitudes smaller compared 
to the original data. A more complete description 
of the developed algorithm and its numerical val- 
idation on high-Reynolds number magnetohydro 



dynam ic turbulence can be found in (jUritskv et al 
2010bh . 



Both data sets (MDI and EUVI) are charac- 
terized by broad distributions of pixel values (Fig. 
[3]) , encompassing more than three orders of magni- 
tude for the magnetic field and about two orders 
of magnitude for the EUVI flux. In an attempt 
to make the event detection consistent with the 
dynamic range of each data sets, we set the detec- 
tion thresholds to percentile levels p = 95.0, 97.0, 
99.0, 99.5 % defining the percentage of image pix- 
els above a threshold. The dimensional thresholds 
corresponding to the chosen set of percentiles in 
MDI and EUVI data are shown with vertical lines 
in the Fig. |3] Table Q] shows the thresholds values 
and the amounts of detected events corresponding 
to these thresholds. 

By applying the described meth od, we identi- 
ficd all distinct path-connected sets (jSherbrooke et al 



19961) of spatiotemporal coordinates satisfying 
the threshold condition and representing separate 
events. Fig. |4] provides an illustration of some 
of the events detected in the conjugate MDI and 
EUVI data sets at p = 95%. 

Each solar event was characterized by the life- 
time T (time interval between the last and the 
first image frames involving the event), the peak 
instantaneous area A covered by the event at any 
given time, the total area S of the entire spa- 
tial region involved in the event during T, the di- 
mension Li a t (Li on ) of this region in latitudinal 
(longitudinal) direction used to estimate the lin- 
ear scale L = {Lf at + L 2 ^) 1 ^ 2 of the event, the 3- 
dimcnsional spatiotemporal volume V of the event 
obtained by summing all district combinations of 
spatial and temporal coordinates attributed to the 
event, and the integrated intensity E of the event 
evaluated by integrating the measured parameter 
(cither LOS magnetic flux or EUVI emission flux) 
over V. 

To quantify multiscale behavior of the detected 
events, the following scaling model has been ap- 
plied: 



/(X)cxX- 
X oc L D " 



(1) 
(2) 



Here, X e {T, L, A, S, V, E} is one of the esti- 
mated event parameters, f(X) is its probability 
distribution, t x is the so-called distribution expo- 
nent of X describing the log-log decay of f(X), 
and D x is the geometric exponent characteriz- 
ing the power-law scaling of the studied parame- 
ter with respect to the linear size L in terms of 
the conditional exp ectation value ~ L Dx 

( Biham et al.ll2001 ). The scaling ansatz described 
by eq. CO)-© is widely used in the theory of 
nonequilibrium critical phenomena, in particu- 
lar, f or describing proc esses of critical percola- 
tion (Munoz et al.l 1999), fractal surface growth 
(jBarabasi fc Stanley! Il995l), an d self-organized 



criticalitv (|Bak et al.lll987L Il988t ) 



Throughout this paper, we use the "hat" sym- 
bol to distinguish coronal scaling exponents D x 
and t x from their photospheric counterparts D x 
and t x . 

Fig. [5] shows the statistical dependence of mag- 
netic and emission event parameters (lifetime T, 
peak area A, total area 5, spatiotemporal volume 
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V and integrated flux E) on the linear size L of the 
event obtained for the 95% thresholds. Ovcrplot- 
ted with the scattcrplots are averaged regression 
dependences (solid lines with error bars) computed 
by coarse-graining the data using exponentially in- 
creasing L-bins. 

The scatterplots in Fig. [5] show a power-law 
increase of the studied parameters with L. The 
power law dependence is distorted at the smallest 
L scales due to the insufficient spatial resolution of 
the data, as well as at the largest scales which are 
typically affected by finite-size effects and/or poor 
statistics (|Robinsonllf99llUritskv et al.ll2006h . To 
obtain robust statistical measures of multiscale be- 
havior, it is important to focus on the intermediate 
range of L scales which is free of such distortions, 
and to make the ranges of all other parameters 
consistent with the selected L range. Based on 
the shape of A(L) dependence, we set the smallest 
and the largest linear scales to respectively Li=3 
Mm and L 2 =12 Mm, with L x < L < L 2 . We 
found that the upper scale limit L 2 can be in- 
creased up to at least f5 Mm without affecting 
the main results of this study. Reducing the lower 
limit L\ leads to less accurate estimates caused by 
the inclusion of small-scale events which are not 
properly resolved. The same L-scale limits were 
applied to MDI and EUVI data sets to ensure con- 
sistent statistical estimates. 

The geometric exponents D x and D x describ- 
ing respectively photospheric and coronal events 
were determined based on power-law regression 
fits to the original (non coarse-grained) data in- 
side the interval L G (Li,L 2 ). The tilted dashed 
lines in each of the panels of Fig. [5] show the log- 
log slopes of the regression fits. The slopes are 
supplemented by the average values and standard 
errors of the evaluated geometric exponents. The 
obtained regression models were used to calculate 
scaling ranges of all other parameters; these ranges 
were subsequently used to measure the distribu- 
tion exponents: 

l£ (a x Lf m , a x L®^ for the photosphere (3) 
X e (a x Lf x ,a x L2 m ) for the corona, (4) 



where a x and a x are the regression coefficients 
evaluated for the two solar regions and L\ (L 2 ) are 
the lower (upper) limits of the studied L range as 
explained above. The studied ranges of L scales 



are shown in each of the Fig. [5] panels by verti- 
cal dot-dashed lines, with the matching intervals 
of scales of conjugate parameters computed based 
on eq. © shown by horizontal lines. 

Probability distributions of MDI and EUVI 
events for all four percentile thresholds are pre- 
sented in Fig. [6l The ranges of scales used to 
evaluate power law slopes of each histogram were 
calculated using eq. [3] These ranges are marked 
by vertical dash-dotted lines and correspond to 
the horizontal lines in the previous figure. It can 
be seen that the distributions decay as approxi- 
mate power laws, consistent with earlier observa- 
tions of po wer law statistics in solar photosphere 
and corona dCrosbv et al] 1993 j_ Abramenko et al 



20021 ; lAschwanden fc Parnelll l2002t lAschwanden 



201ll ICrosbvll201lt lAschwandenll2f)13l) 

The broad-band power-law behavior of MDI 
distributions suggests that the dynamics of the 
LOS photospheric magnetic field is an inhe rently 
multiscale process (jUritskv fc Davilall2012h . Our 
threshold-based definition of magnetic events is 
clearly not limited to the injection of new magnetic 
bipoles into the photosphere and their subsequent 
submergence into the convection zone. In essence, 
we deal with an apparent dynamics associated 
with crossings of the percentile detection threshold 
by the local magnetic field. The resulting sets of 
events likely contain a combination of truly emerg- 
ing and submerging magnetic bipoles as well as 
magnetic flux coalescence, fragmentation, cancel- 
lation, injectio n of upwardly propagat ing Alfvcnic 



{)agatn 
20071 ). and other 
effects leading to spatiotemporal variability of the 
photospheric field. As we demonstrate below, 
the dominating factor leading to the multiscale 
statistics of this complex set of events is random 
footpoint shuf fling driven by tur b ulent convective 



flows (s ee e.g. lLamb et al.l (120081) :IChaouche et al 
(l20Ill) : ICrouch et al.1 (j2007h : ISimon et al.1 <|l995L 
200lh ). 



The values of EUVI distribution exponents, as 
well as SO HO EIT exponents r eported in our ear- 
lier study (jUritskv et al.1 120071 ). suggest that so- 
lar corona operates in a highly intermittent turbu- 
lent state consistent with models of self-organize d 
criticality (|Bak et al.llf987tlLu fc Hamiltonlll99ll) . 
It should be noted, however, that the exponents 
reported here are obtained using spatiotemporal 
definition of events which is conceptually much 
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closer to measuring avalanches in numerical simu- 
lations of self-organized criticality than most of 
the definitions used in previous works on flare 
statistics conducted by other authors - except 
for a few case studies faerghmans et al. 1998 ; 
Berghmans fc Clettd I1999T focusing on specific 
coronal conditions. The energy distribution ex- 
ponent te < 2 indicates that the coronal dynam- 
ics is likely to be dominated by large events as 
opposed to P arkers scenario of nanoflare heating 
()Parkerlll988h . Coronal energy scaling exponents 
below the value 2 have been obtained in many pre- 
vious studies including those based on spatiotem- 
poral detection of coronal brightcnings and their 
spatial detection with subsequent integration of 
the emission fluxes over fixed time interval (e.g 



Aschwanden et al. 201 )U a In. As we argue below, 



such exponent values suggest that small heating- 
events are not frequent enough to account for the 
bulk coronal heating. The difficulty with ascrib- 
ing this role to the large events is that they arc 
unable to supply the require d energy loss rates of 
10 5 - 10 7 crg/cm 2 /s (see e.g. lAthavl[l966h . 

The exponent tt is consistent with previous 
analyses of threshold-dependent inter-occurrence 
times of x-ray burs ts measured over the whole Sun 
( Baiesi et al.ll2006r ). as would be expected if total 



emission were a sum ove r individual avalanches 
( Corral fc Paczuski 1999h . We also note that 
this and other EUVI distribution exponents are 
in an approximate agreement with the results of 
forc ed MHD simulation of tur bulent coronal heat- 
ing (jDmitruk fc Gomezl[l997h . 

Table [2] provides a summary of geometric and 
distribution scaling exponents in the photosphere 
and corona for all four detection thresholds. In 
addition, the table contains the auxiliary scaling 
indices A x = D x (t x — 1) and A^ = D x (t x — 1). 
For a given threshold and data type, these indices 
are expected to be independent of X (constant 
across all event parameters) due to the probabil- 
ity conservation implying f{X)dX ~ f(L)dL for 
any X. A deviation of (A^ ) from the values 
D l {t l -1)=tl-1 [Dl{tl - 1) = t l - 1) by 
more than ~ 10% signals an insufficient statistics 
and/or an inadequate scaling model. Such cases 
are marked by asterisk in Table [5] 

Table [2] shows that the probability distributions 
of MDI and EUVI events corresponding to the 
same percentile level are characterized by signif- 



icantly different power law slopes yielding distinct 
t x exponents. In particular, the t x values charac- 
terizing photospheric events obtained at the 95% 
detection threshold are substantially higher com- 
pared to the corresponding t x values describing 
EUVI events obtained at the same threshold. 

We found that the difference in the scaling be- 
havior of the two solar regions originates in the 
difference of the intermittency level of events de- 
tected at the same percentile thresholds. As Ta- 
ble Q] indicates, for any given percentile thresh- 
old, the number of MDI events is considerably 
larger compared to the number of EUVI events. 
For the lowest applied threshold, the discrepancy 
reaches an order of magnitude. This suggests that 
for the same p, the magnetic events are much 
more intermittent than the coronal events. Ap- 
plying a smaller threshold to a spiky MDI magne- 
togram results in much more threshold crossings, 
and hence more time-integrated events. Due to 
a much smoother spatial pattern of the coronal 
emission field, the same procedure leads to a less 
significant increase in the number of the detected 
coronal events. 

Physically, it is clear that comparing popula- 
tions of MDI and EUVI events obtained at the 
same detection threshold is not justified because 
such events represent different types of conditions 
characterized by incomparable occurrence rates. 
Indeed, there is no reason to expect that a fre- 
quently occurring photospheric condition is de- 
scribed by the same statistics as a rarely observed 
class of coronal events. To make the compari- 
son of the photospheric and coronal events more 
consistent, one needs to make sure that the de- 
tected populations of events have approximately 
the same size - in other words, they are described 
by the same degree of intermittency. To attain this 
goal, one needs to set the magnetic threshold to 
a higher percentile level compared to the emission 
threshold. 

Based on the population sizes provided in Ta- 
ble [U the 99.0% (99.5%) magnetic events have 
approximately the same degree of intermittency 
as the 95.0% (97.0%) heating events. Pairwise 
comparison of these sets of events (99.0% MDI 
vs. 95.0% EUVI, 99.5% vs. 97.0%) shows that 
they are described by drastically different geomet- 
ric exponents but surprisingly close distribution 
exponents, especially those for the time-integrated 



G 



parameters S, V and E (see Table [2]). As we 
show in Section 13.31 the distinct values of the lin- 
ear scale distribution exponent in the photosphere 
and corona are consistent with the geometries of 
the intermittent events in the two system. 

Fig. [7] presents a direct comparison of the 
probability distributions of MDI and EUVI events 
over the three measures S, V and E involving 
spatiotemporal integration. The distributions are 
built using the pairs of percentile thresholds yield- 
ing consistent population sizes. To construct the 
distributions, we consider only the events whose 
linear sizes lie within the scaling interval of inter- 
est (Xi, L2). No other transformations besides the 
standard probability density normalization were 
applied to the histograms. The obtained MDI 
and EUVI distributions arc remarkably similar 
in shape and are described by fairly close prob- 
ability values. Within the statistical accuracy, 
they effectively coincide for both combinations of 
thresholds. The plotted spatiotemporal measures 
of coronal emission events obey the same occur- 
rence probabilities as the corresponding measures 
of photosphcric magnetic events. 

The similarity of coronal and photospheric dis- 
tribution exponents is also evident from Fig. [8] 
which shows time evolution of the exponents eval- 
uated within a 90-minute moving window. The 
wide error bars reflecting the lack of the suffi- 
cient number of events prevent us from drawing 
solid conclusions regarding the observed trends, 
although it seems plausible that the probabilistic 
laws describing the two data sets are dynamically 
coupled. 

These observations are in line with the analy- 
sis of unstable magnetic discontinuities identified 
based on linear force-free extrapolation o f pho- 
tospheric field dVlahos fc Georgoulisl 120041 ) Us- 
ing this technique, Vlahos fc Georgoulisl ( 20o3 ) 
have found a large number of magnetic discon- 
tinuities whose free magnetic energies and vol- 
umes obey power-law distribution functions. The 
log-log indices of the distributions were found to 
be in agreement with the indices describing the 
occurrence frequencies of solar flare energies re- 
port ed in the literature. Based on this compari- 
son, Vlahos fe Georgoulisl ( 2004 ) speculated that 
the statistics of flares may result from a preexist- 
ing free energy fragmentation in the solar corona. 
Solar flares occur in unstable coronal volumes as- 



sociated with discontinuities in the coronal mag- 
netic field. It is expected that not every single un- 
stable magnetic volume generates a flaring event, 
but most of the flares do occur in such pre-existing 
unstable volumes containing strong nonpotential 
magnetic field. 

Our statistical analysis based on conjugate sets 
of coronal and photospheric observations provides 
an opportunity for a more direct test for this cou- 
pling scenario. Taken at face value, the statis- 
tics shown in Fig. [7] support the picture in which 
the coronal emission events are generated in a 
random subset of unstable magnetic volumes cre- 
ated by the photosphcric motion. The quantita- 
tive analysis of scaling relations between coronal 
and photosphcric power-law exponents presented 
in Section 13.31 below lead us to the same conclu- 
sion. At the same time, our results do not rule out 
a possibility of a different underlying physics in- 
volving multiscal e pulses of upwardly propagating 
Alfve nic waves dSuzuki fe Inutsukal l2006t lUlrichl 



19961 : iFuiimura fc Tsunetal 12001 as discussed in 
section 13.41 



3. Discussion 

In this section, we look closer at the spatiotem- 
poral geometry of magnetic and emission events 
in an attempt to derive scaling relations between 
the observed scaling exponents. The analysis is 
validated using the averaged scaling exponents de- 
scribing MDI events detected at p = 99% and p = 
99.5% thresholds, as compared to the averaged ex- 
ponents of EUVI events obtained at p = 95% and 
p = 97%. As mentioned earlier in the text, the 
sets of events corresponding to these thresholds 
are described by comparable levels of intermit- 
tency resulting in consistent population sizes. The 
exponent values used for the validation and are 
provided in Table |3l While the numerical values 
of the exponents vary with the detection thresh- 
old (Tabled]), the relations between the exponents 
discussed in this section are expected to be insen- 
sitive to the threshold as long as the MDI and 
EUVI events have comparable population sizes. 
As before, we use the hat symbol (") to distin- 
guish coronal parameters from their photospheric 
counterparts. 
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3.1. Scaling of the magnetic events 

According to the accepted scaling model ((2|), 
the scaling of linear size of the event with the event 
duration is governed by the geometric exponent 
D T : 

L oc T^ Dt = T H , (5) 

Since D T w 2.18 (Tabic E]) the diffusive expo- 
nent H w 0.46 is close to the value 1/2. This 
means that the linear size of the perturbation 
scales approximately as a square root of the 
elapsed time. Such scaling is characteristic of 
the cl assical (as opposed to fractio nal) random 
walk ([Mandelbrot fc VanNessl Il968h which can 
be a manifestation of random footp oint shifting 
of the phot o spheri c field (see e.g. Simon et al.l 
(|l995l . l200lh : iRastl (l2003h : ICrouch et al.1 foQQm . 
The idea of the random footpoint shuffling has 
been recently used for constructing a statistical 
fractal-diffus ive self-organized c riticality model 
(FD-SOC) (lAschwandenl l2012allbh The stan- 
dard FD-SOC model |Aschwanden| l2012bh with 
H = 1/2 predicts the geometric and distribution 
exponents which are fairly close to the exponents 
of the MDI events reported in Table [3J 

Assuming that the average magnitude of the 
horizontal convective flow associated with the 
shuffling is independent of L within the studied 
range of spatial scales, wc get a simplified picture 
(Fig. E{a)) in which the swept area of the event 
scales approximately as S oc £ T, where £ ~ A 1 / 2 is 
a measure of the instantaneous linear scale_of the 
event. Therefore we have 5" oc A X I 2 T oc L~ VDt 
which, taking into account S oc L Ds , yields 



D, 



Da 
2 



Dn 



(6) 



By plugging in the measured values of Da and 
Dt we get Ds ~ 2.88. The direct evaluation of 
D$ based on the S(L) regression leads to Ds = 
2. 78 ±0.11 showing that the relation © holds with 
a reasonable statistical accuracy. 

The exponent Dy describing the scaling of spa- 
tiotemporal volume V of magnetic events can be 
also derived in frames of the random shuffling 
model implying that V ~ A T (since the cross- 
section area A is orthogonal to the time axis, see 
Fig. Ha)). Making use of A oc L Da and T oc T Dt 
we get 

Dy = Da + Dt- (7) 



The relation (|Tj) predicts Dy 3.57 which agrees 
with the observed value Dy = 3.43 ± 0.12. 

It is important to clarify that the equation (|6|) is 
derived assuming that the instantaneous shape of 
magnetic elements is purely 2-dimensional. While 
this assumption is certainly not justified for high- 
resolution magn etograms featuring ir regular mag- 
netic geometry (jCadavid et al it seems to 
be acceptable for describing SOHO MDI struc- 
tures covering up to 10 Mm 2 at a time (the largest 
A consistent with the chosen L range, see Fig. [5]) 
since such structures are too small to display iden- 
tifiable fractal boundaries at the available resolu- 
tion. 

The dependence of the instantaneous linear 
size £ of magnetic elements on lifetime T reflects 
the dynamic behavior of the underlying physi- 
cal processes acting at small spatial scales, such 
as e.g. magnetic flux coalescence, fragmenta 



tion, emergence and subm e rgence (|Lamb et al 



20081; IChaouche et al.l lioTl ICrouch et al.l l2007t 



Simon et al.1 Il995l l200li IWelschl l2006h . Even 



though these processes are not fully resolved in 
the studied data, their effective contribution can 
be evaluated based on the temporal scaling of mag- 
netic element sizes as expressed by the measured 
geometric exponents. By combining £ oc A 1 / 2 , 
A oc L Da , and L oc T 1 / Dt we obtain 



£ocT Hl , H l 



Da 
2D t 



(8) 



Here, the diffusion index Hi characterizes a small- 
scale transport of the photospheric magnetic 
field. Substituting D A m 1.39 and D T w 2.18 
we get Hi w 0.32. The obtained index shows 
that the small-scale magnetic field convection is 
strongly subdif fusive (jBouchaud fc Georges 19901 : 
Isichenkolll992 ) , meaning that the size of the mag- 
netic elements varies with time significantly slower 
than it would be expected for a regular diffusion 
described by a classical Brownian walk defined 
by Hi — 1/2 (linear displacement proportional to 
the square root of the elapsed time). A subdiffu- 
sive transport usually takes place when the stud- 
ied system contains spatial subdomains where the 
transported quantity bec omes partially "trapped" 
(jMetzler fc Klafterl |2000h ■ In solar photosphere, 
such trapping can be due to the small-scale tur- 
bulence accompanyin g the granulation process 
(jEspagnet et al.lll993D which is closely related to 
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the s urface dynamo problem (jVogler fc Schiissler 
l2007h . According to Kolmogorov's refined sim- 
ilarity hypothesis, the third moment of the dis- 
tribution function of correlated velocity incre- 
ments Au observed at the characteristic dis- 
tance £ scales as £ 3Hl (jStolovitzky fc Sreenivasan 



19941 ). Comparing this expression with Kol- 



mogorov ' s exact relation for inertial range tur- 
bulence ( Monin fe Yagloml 1197.4 iBiskamnl l20f)l 



(Au(£) 3 ) = — I (e) I (e being the dissipation rate 
per unit mass) one gets Hi = 1/3. This pre- 
diction is fairly close to the Hi -value obtained 
for the photosphere in our study showing that 
that the range of sizes of magnetic elements ad- 
dressed by this analysis (£ = A 1 ! 2 ^2 — 4 Mm, see 
Fig. [5] ) is likely to fall into the i nertial range of 



the g ranulation-driven turbulence (jEspagnet et al 



119931 ). and if so, the photospheric turbulence is lo- 
cally isotropic. 

3.2. Scaling of the emission events 

In contrast to the photospheric events evolving 
approximately as L oc T 1 / 2 , temporal dependence 
of linear sizes of coronal events is strongly su- 
perdiffusive as revealed by H = 1/Dt ~ 0.74. 
The small-scale coronal diffusive exponent Hi = 
Da/(2Dt) ~ 0.59 is also larger than the classical 
diffusion value. Superdiffusive behavior is known 
to arise in strongly nonlinear media containing 
spatially distributed regions in which the diffusive 
trans port is anomalously fast (jMetzler fc Klafter 
2000). Compared to the classical random walk 



with temporally uncorrelated increments, the su- 
perdiffusive walk, such as the one modeled by frac- 
tiona l Brownian motion ([Mandelbrot fc VanNess 
19681 ). features significant positive correlations 
across all time scales involv ed in the process 
(|Turcottd Il989t iFederl Il988h . The fact that 
H ^ Hi suggests that the dynamics of coronal 
emission regions is a sum of two types of motion 
- a superdiffusive expansion of the event area and 
a displacement of the center of the event over a 
comparable distance (Fig. EJb)). 

The expansive superdiffusive growth of emis- 
sion regions could be a manifestation of energy 
avalanches develop i ng in the co r ona (see e.g. 



potential magnetic field configurations existing 
in the corona release their free energy through 
a chain interaction of multiple spatially local- 
ized instabilities s uch as those associated with 
nanofiare heating (IParkerl 1198 1 iKlimchukl [20061: 
Viall fc Klimchukll201ll) . " The perturbation pro- 



ceeds explosively by involving a growing number 
of unstable plasma regions, similarly to the inter- 
action of rollin g grains of sand leadin g to a sand 
pile avalanche dBak et al.lll987l Il988l ). This sce- 
nario, known as self-organized criticality (SOC), 
has been extensively studied in the context of 
so lar flares, starting fro m the pioneering work 
of iLu fc Hamilton! (|199lf ) and followed by a va- 
riety of more recent num erical simulations (see 



Charbonneau et al.l (|2001l ) for a brief review) . It is 



important to note that coronal energy avalanches 
may or may not be topologically compact. In the 
latter case, the super-diffusivity of solar corona 
can be a manifestation of multiple dissipation re- 
gions a ppearing in respon se to the same initial 
trigger (jTorok et al.l 1201 lh . The ensemble-based 
behavior of these regions would involve spatial 
jumps (such as those modeled by Levy flights) 
associated with the appearance of new centers of 
activity, leading naturally to a superdiffusive H 
value. The evolution of individual regions could 
follow a different scaling law and be subdiffusive 
during the impulsive phase of flares ( Aschwandenl 



2012al ). 



Despite their major success in explaining scale- 
free cooperative phenomena, most of the SOC 
models are not solvable analytically (|jensenlll998l ) 
making it difficult to obtain closed-form rela- 
tions characterizing observed statistical scaling 
laws. The first-principle analysis of such mod- 
els is often an intractable task due to the in- 
volvement of higher-order statistical moments 
(jVespignani fc Zapperil Il998l : ISethna et ail l200ll ) 



avalancnes develop i ng m tne co r ona (see e.g. 

Crosby et al.l (ll993l);ICrosbyl (l201ll ) ;ICharbonneau et alj 



(2002) 



2007D ) 



( 2001 ); Aschwanden et al. ( 2000bl):lAschwanden fc Parne 



: Morales fc Charbonneau! (|2008h : lUritskv etakl 
It has been suggested that the non- 



requiring full dynamic renormalization group 
treatment (Chang 19921 ). However, this theoreti- 
cal complication could be of little importance for 
the solar activity since its dynamics is not a pure 
SOC but combines both driven and self-organized 
components. 

If the behavior of solar corona was fully analo- 
gous to SOC models, the spatial distribution of 
I^jfree magnetic energy in the corona would have 
to be solely formed by previous flaring events. 
However, coronal dynamics is fundamentally dif- 
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fcrcnt. Although local triggers of coronal instabil 
ities are largely unpredictable ( Aschwandenl 2011 



yields 



20131 ). spatial distribution of free magnetic en- 
ergy supplied by the photosphere is far from ran- 
dom. For this reason, solar corona is a driven 
rather than a fully self-organized avalanching sys- 
tem. The present study confirms this interpre- 
tation by showing that solar photosphere places 
constraints on the occurrence probability of coro- 
nal events of different sizes. 

As the next step, we look further into the 
growth and decay dynamics of coronal emission 
events as characterized by geometric and distri- 
bution exponents. Avalanching dynamics can be 
seen as a special case of spreading dynamics in 
active nonlinear media wit h spatially distributed 
energy sources and sink s ( Vespignani fc Zapperil 



1998 : Munoz et al. 19991 ). Spreading experiments 
provide the most accurate determination of the 
critical poi nt of systems with mult i ple a bsorb- 
ing states ( Grassberger fe Delatorrd 1979() such 
as a valanching systems at self-organized critical- 



ly ( Vespignani et al. 2000h . A small perturba- 



tion representing an initial instability is created 
at the origin of an otherwise absorbing configu- 
ration, leading to a spread of activity. At the 
critical point associate d with power-law distribu- 
tions of avalanche sizes ( Bak et al~lll987 ; Robinsonl 
1994 : Sethna et al. 200lb . the mean-squared de- 
viation from the origin R 2 , the survival prob- 
ability P, and the number N of unstable sites 
scale with time t as R 2 (t) oc t z , P(t) oc t~ s , 
and N(t) oc where z, 5 and f] are the spread- 
ing exponents describing the growth and decay 
of energy avalanches. In our measurements, P 
has the meaning of the cumulative distribution of 
avalanche lifetimes: P(t) = f* f(T)dT oc T ~ fT+1 
and hence 5 = tt — 1 (|Uritskv et al.l 12007 : 
Morales fc Charbonneaull2008l ). 

The number of individual unstable e vents (top 



pling events in sandpile simulations (|Bak et al 
19871 ). nanofiares in the Park er's heating see 



nario (jCharbonneau et"aL 200lh ) associated with 
an aval anche with a lifet i me T scales as oc 
T i+fi+s dMunoz et all Il999t lUritskv et al.l l2007t 
Morales fc Charbonneau 20081 ). Approximating 
this number by the integrated emission flux E we 
obtain E oc T 1+f i +& = zA(i-H}+5) K jfin which 



6= -A 



De 



(9) 



Munoz et al. I (Il999h have argued that the rela- 



tions between the spreading and the distribution 
exponents depend upon the universality class of 
the avalanching system. Here we accept their in- 
terpretation in terms of the directed and dynami- 
cal percolation classes govern ed by the well - known 



Munoz et all (|l999h 



relation zd/2 = r) + 26 (see 
and references therein) which in our case can be 
rewritten as 



zd/2 = ^ 
D T 



tt 



(10) 



where d is the Euclidean dimension of coronal 
avalanches. Also, from the definition of the z ex- 
ponent it follows that 



D T 



(11) 



where we took into account that R 2 ~ S and 
t ~ T oc L Dt . By solving eqs. ([TOllTTj) for d 
and substituting the measured coronal exponent 
values we find that d ~ 3, with the value z ks 1.47 
predicted by cq. (|ll[) being in an agreement with 
the beha vior of three-dimen sional dynamical per- 
colation (|Munoz et al.lll999l ). 

These estimates suggest that coronal dissipa- 
tion takes place in 3-dimensional plasma volumes. 
Elementary instabilities associated with individ- 
ual nanofiares can, in principle, have a different 
dimensionality, for instance d = 2 for reconnect- 
ing current sheets. However, our estimates show 
that the interaction network controlling growth 
and decay of coronal emission regions is likely to 
be embedded in a 3-dimensiona l space, and is also 



space-filling (|Mandelbrotlll977l) 



is often taken for granted ( 


Aschwanden et al. 


2000b: Parnell & Jupp 20001 


iKrucker & Benz 


1998: Mitra-Kraev & Benzl 200lh. The analv- 



sis conducted here justifies this choice speak- 
ing in favor of a three-dimensional emitting 
plasma volume, as opposed to more complex ge- 
ometries described by d < 3 as predicted by 
some statistical-physical model s of solar flares 



(jMelntosh fc Charbonneaull200lD . 
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The derived three-dimensional geometry of en- 
ergy avalanches in the corona is consistent with 
other scaling exponents. The swept area S of each 
event can be considered as a two-dimensional pro- 
jection of a d-dimensional shape. The dimension 
dp of the projection (approximated in our case by 
Ds by the de finition of fractal dimension (e.g., 
Turcottd[l997h ) is related to the dimension of the 
projected object as d p = d — 1 (jMandelbrotil 1 982[) . 
For d = 3, the relation predicts Ds = 2 which 
matches exactly the measured value 2.00 ± 0.03 of 
this coronal exponent. Also, assuming (as before) 
that the integrated emission flux E is proportional 
to the number of the unstable coronal sites con- 
stituting an avalanche, we expect E oc L De with 
De = 3, which is in an agreement the empirical 
result D E = 3.06 ± 0.05. 

3.3. Stochastic coupling relations 

Our observations indicate that the MDI and 
EUVI events are described by rather similar prob- 
ability distributions of S, V and E parameters 
but, at the same time, by considerably different 
geometric scaling laws yielding distinct sets of D x 
and D x exponents (see Tablc[3]). This section aims 
to resolve the apparent controversy by connect- 
ing the two groups of exponents in through a set 
of stochastic scaling relations between the photo- 
sphere and corona. In section these relations 
are used to provide a new insight into the physics 
of the underlying photosphere - corona coupling. 

As a starting point, consider the subsets of coro- 
nal and photospheric events characterized by the 
same probability of occurrence. Since the prob- 
ability distributions of the swept areas of coronal 
and photospheric events coincide, this requirement 
is equivalent to S — S (coronal measures are la- 
beled by the hat symbol). 

Eq. © suggests that the linear scales describ- 
ing the considered events are related as L Ds oc 
L Ds yielding 



L oc L a , 



D s 



1.39. 



(12) 



meaning that the apparent linear sizes of the coro- 
nal events are larger for the same swept area. The 
new scaling index a defined by the ratio of the 
measured values of Ds and Ds exponents couples 
the linear scales of equally probable coronal and 



photospheric events. The fact that a > 1 means 
that the characteristic length of coronal emission 
regions grows faster with the swept area compared 
to the length scale of photospheric events. The 
difference in the scaling behavior of the two sys- 
tems expressed by a should hold under any linear 
transformation maintaining S oc S. The power- 
law relation (fT2]) between probabilistically equiv- 
alent photospheric and coronal scales accounts for 
nonlinear corrections to such geometric transfor- 
mations associated with a physical interaction be- 
tween the two system. 

Next, consider the lifetimes of the selected sub- 
sets of events described by S — S. Since T oc L Dt 
and T oc L Dt , we have a similar scaling law for the 
lifetimes of the coronal and photospheric events, 
with a differing value of the power-law index: 



f ocT' 3 , P 



DsDt 
D T D S 



0.87. 



(13) 



The obtained value of temporal coupling index 
/? indicates that the duration of coronal emission 
events increases slower with the event area com- 
pared to the duration of photospheric events. 

The relations (fT^]) - (fT^j) suggest that the ob- 
served similarity of f(S) and f(S) distributions 
is achieved through a nontrivial rescaling of spa- 
tiotcmporal dimensions of solar events, with the 
coronal emission events lasting a shorter time 
but propagating over larger distance compared to 
the photospheric shuffling events described by the 
same occurrence rate. 

Our goal now is to derive coupling relations ex- 
pressing coronal geometric and distribution expo- 
nents as functions of the photospheric exponents 
and the empirical indices a and /?. 

By definition, the linear scale geometric expo- 
nent is identical to unity in both data sets, and 
so 

D L = D L (=1). (14) 

The exponent Ds is already connected with Ds 
through the empirical scaling index a (|12[) . It can 
be also linked to other photospheric indices. By 
combining the swept area expression for the pho- 
tosphere © with the coupling relation (fT2")l we 
obtain S oc (i" Q )^ L+Dr . After equating S with 
S oc L Ds , this equation provides an alternative 
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form for the coronal exponent D$'- 



D, 



Da + 2£> T 
2a 



(15) 



This relation predicts Dg ~ 2.0 which coincides 
with the measured value for this exponent (2.00 ± 
0.03). 

Similarly, by plugging (IT2|) and (|T3l) into the 
scaling ansatz T oc L Dt and substituting L Dt for 
T we obtain 

£> T = -D T (16) 



which yields the correct result Dp ~ 1-36 as ex- 
pected from the definition of a and /3. 

For deriving the exponent Da describing the ge- 
ometric scaling of the peak area A of coronal emis- 
sion events one should take into account fractal 
(non- integer) dimension cLf of the emission regions 
defined by A oc £ dF . In frames of our model (Fig. 
Kb)) the S covered by an emission event re- 
sults from the combination of the avalanche-like 
expansion of the unstable region and the system- 
atic displacement of its ccntroid. The charac- 
teristic linear scales of these effects are given re- 
spectively by I and L from which it follows that 

S oc LI = L d F and therefore 



D A 



Dc 



1 



(17) 



The obtained relation allows us to evaluate the 
fractal dimension dp « 1.6 without performing 
direct scaling analysis which would likely to be in- 
accurate due to the smallness of i and limited spa- 
tial resolution of EUVI images. Using (JTTJ) along 
with the coupling relation for D$ (jT5"|) we get 



( — -1 

a 



(18) 



By the definition of a and dp, the obtained expres- 
sion matches the measured value Da = 1.60±0.02. 

The exponent Dy can be related to Dy by 
making use of the similarity of probability dis- 
tributions of spatiotcmporal volumes of coronal 
and photospheric events (Fig. [7]), by analogy with 
the analysis of the exponent Dg presented above. 
For a given occurrence rate, V = V and hence 
L v oc L Dv , which after applying the spatial 



rescaling (IT21) provides 



Dy 

a 



(19) 



Comparing (|15j) and (fl9|) leads to the identity 
(DsDy)/ (DyDs) = 1 which is approximately 
fulfilled (substitution of the measured exponents 
yields -0.95). 

Expressing the emission flux exponent De via 
the photospheric exponents presents a critically 
important (albeit somewhat more challenging) 
task because this exponent, along with the dis- 
tribution exponent te, can be used for forecasting 
coronal emission level. Using the relations (|TU)) 
and (fTTj) for the coronal spreading exponent z 
along with the coupling relations (fT5")) and (|16l) . 
after some algebra we arrive at 



Di 



$D T ( d D & 



\2pD 7 



(20) 



Note that this is an intermediate result because it 
depends on the coronal exponent vr which is yet 
to be related to derived. 

The general form for any distribution exponent 
f x follows directly from the probability conser- 
vation condition f(X)dX — f(X)dX where X 
and X stand for the same type of measures in 
the corona and photosphere. Assuming a power- 
law scaling relation X = X 7x between the coro- 
nal and photospheric measures and applying the 
power-law scaling model ([T]) we get f(X)dX oc 
X—i^'Xt'-HX oc X- T *dX and consequently 



1 D x 

- + 1, 7x = -^a (21) 



7x 



where the expression for j x follows from X oc 

lP* - (L a ) 3 - - ^x l ! D *) a ) b * - X^ a . By 
using the expressions for coronal geometric expo- 
nents derived above, one can easily calculate the 
corresponding value of j x and on this basis f x . In 
particular, for the linear size distribution exponent 
we have jp = a (since Dp = Dp = 1) and 

TP - 1 



TP 



1. 



(22) 



In its turn, the lifetime distribution exponent can 
be derived in a similar fashion by taking into con- 
sideration eq. tfTB"]) which yields 7t = /? and there- 
fore 

f T = ^_i + L (23) 
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The relations ([22)1 and (|23|) result in the estimates 
fx ~ 2.55 and tt ~ 1-91 which fall within sta- 
tistical uncertainty of the measured values tl = 
2.57 ± 0.13 and f T = 2.02 ± 0.07. 

Finally, after the substitution of (J23J) into ((20)) 
we obtain the geometric exponent of the emission 
flux as a function of the photosphcric scaling ex- 
ponents and the empirical coupling indices a and 



D, 



dD s /2 + {(3 - t t + 1)D T 



(24) 



For space-filling 3-dimensional emission regions 
(d = 3) the relation jH} predicts D E ~ 3.12 
which is approximately consistent with the ob- 
served D E = 3.06 ±0.03. 

The last column in Table [3] provides a summary 
of the coronal distribution exponents predicted us- 
ing the scaling relations discussed in this section. 

3.4. Energy considerations 

The analysis of scaling relations presented so 
far reconciles the obtained results showing that 
the difference in the geometry of photospheric and 
coronal events is expected for physical reasons, 
and is consistent with the similarity of the statis- 
tics describing spatiotemporal event measures S, 
V and E. The derived scaling laws suggest a 
simple coupling scenario in which the footpoints 
of the coronal flux tubes undergo a quasi one- 
dimensional horizontal convective motion given by 
a vector sum of uncorrelated displacements (a ran- 
dom walk behavior) while the coronal emission 
events are three-dimensional, with the grow and 
decay dynamics being similar to that of SOC mod- 
els. When integrated over space and time, the two 
groups of events exhibit indistinguishable distribu- 
tion functions. 

The physical agent that provides the most natu- 
ral explanation for the observed probabilistic cou- 
pling between solar photosphere and corona is the 
nonpotential magnetic energy supplied by the pho- 
tospheric shuffling. 

The normal Poynting flux S n of energy through 
the base of the corona in the ideal magnctohydro- 
dynamic approximation can be written as 



x B 



x B 



(25) 



the unit normal vector (for simplicity the mag- 
netic diffusivity is assumed to be negligible). A 
similar expression can be written for the potential 
magnetic energy flux; the difference between the 
two fluxes yields the flux of f ree magnetic energy 
into the corona (|Welschll200i : 



Q(F) - _ (r, 



B 



(P. 



(v„B h - v h B n ) (26) 



in which v is the photospheric velocity and n is 



where the subscripts h (n) denote the horizontal 

(normal) components of the velocity and magnetic 

(P) 

field vectors and B^ is the potential horizontal 
field. The studied FoV is close to the center of the 
solar disc where SOHO MDI flux counts provide a 
reasonable proxy to B„. The behavior of the hor- 
izontal velocity can be evaluated indirectly based 
on the tracks of the magnetic features as discussed 
in section [3l 

Assuming that the difference between the to- 
tal and the potential horizontal magnetic fields re- 
mains constant during each event and ignoring the 
contribution from the u„B/, term which cannot be 
evaluated from the studied data, we find that on 
average, 

S<p cc (v h ) (B n ) (27) 

with (vh) and (B n ) being the characteristic values 

of the horizontal flow velocity and the line of sight 

(p) 

magnetic field, respectively. The sign of S n in eq. 
(|27p is positive since we are not interested in the 
events producing a downwardly directed Poynting 
flux S { n F) < 0. 

According to our statistical results, the dynam- 
ics of the photospheric footpoint shuffling is close 
to classical random walk with H = 1/2. This 
dynamics has been shown to be critically impor- 
tant for the generation of the upward Poynting 
flux, especially in strongly magnetized intergran- 
ular lanes experienc ing horizontal vortex motions 
( Shelvag et al.112012 ). The velocity field satisfying 
the H = 1/2 condition is described by a stable 
probability distribution w ith finite and c onstant 
first and second moments (|Turcottelll997l ). Based 
on this model, one can assume that (vh) remains 
constant during each event and is not a function 
of its linear scale L (such a dependence would im- 
ply a super- or sub-diffusive behavior of the mag- 
netic events which is not observed). The quan- 
tity (B n ) can be easily obtained based on the con- 
ducted scaling analysis. The spatiotemporal sum 
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of magnetic flux counts during the events scales as 
oc L De while the spatiotemporal volume scales as 
oc L Dv . The ratio between these two quantities 
gives the required average field: 



(B n ) oc L 



D E -D^ 



(28) 



The free magnetic energy W produced by any 
given shuffling event can be evaluated by integrat- 
ing eq. (|2"T)l over the involved surface and the 
event lifetime: 



W OC (V h ) / dt 



dxdy 



(29) 



x,y£a(t) 



in which a(t) is the instantaneous area of the event 
at time t. It is easy to verify that the domain of 
the integration in cq. (|29| is nothing but the swept 
area S of the event, according to our definition of 
this parameter. By making use of this fact and 
substituting ((28]) for the magnetic field we arrive 
at the geometric scaling law for the energy of the 
shuffling event: 



W oc L De - d ^ +Ds . 



(30) 



Substitution of the estimated values of the pho- 
tospheric scaling exponents involved in this rela- 
tion yields W oc L Dlv with D w w 3.15. This 
estimate is in a decent agreement with the geo- 
metric exponent De = 3.04 ± 0.03 characterizing 
the power-law scaling of the coronal emission flux 
with L. Therefore, the shuffling energy and the 
coronal emission energy in the populations of coro- 
nal and photospheric events described by the same 
rate of occurrence are directly proportional. While 
the proportionality itself does not necessarily im- 
ply a causal link, it clearly speaks in favor of a 
photospheric input working as a driver of coronal 
dissipation. 

The described coupling scenario is reminiscent 
of the famous Parker's scenario of nanoflare coro- 
nal heating. However, the similarity can be mis- 
leading for the following reasons. First, the phys- 
ical scales of the coronal events included in our 
analysis is incomparably larger than the size of the 
nanoflares. This is because the detection thresh- 
old used to extract the events was tuned to a fairly 
high emission level which produced small appar- 
ent sizes but in fact ensured that most of the 
low intensity coronal brightenings were excluded 
from the statistics. Secondly, if we postulate that 



the smaller events (unresolved by our technique) 
obey the same statistical laws as the observed 
events, the obtained value te ~ 1.5 < 2 of the 
emission exponent predicts that the energy out- 
puts from the small events (possibly including the 
nanoflares) is negligible compared to the contri- 
bution of large events since the ensemble averaged 
emission flux 



(25) 



Ef(E)dE oc E 



2-T E \E„ 

\e„ 



(Ejnax ) 



0.5 



(31) 

is controlled by the largest flares. On the other 
hand, the ratio between the free magnetic energy 
of a given shuffling event ([30| to the average total 
magnetic energy per unit area during the event 
increases with decreasing L: 



W 



Free magnetic energy 
Total magnetic energy J B\ dv 

L D E -D v + Ds £ 3 - 15 



£2(£>e-£) v )+3 



L 



3.74 



L 



-0.84 



(32) 
(33) 



in which we used the energy scaling relation ([30]) 
for W, equation (|2"5)l for the average LOS mag- 
netic field during the event, and have assumed a 
3-dimcnsional coronal volume. The fact that the 
resulting power law exponent is negative indicates 
that the small-scale magnetic events are on av- 
erage much more effective as free energy sources 
compared to the large events. The higher rela- 
tive occurrence rate of the small shuffling events 
could make this dependence statistically quite im- 
portant. However, the net free energy outcome 
from the small photospheric events is still smaller 
than that from the large events because the power- 
law exponent Tw of the W distribution predicted 
by the probability conservation (Dw{t~w ~ 1) = 
Dl(tl — 1) => Tw ~ 1-7) is, again, below the 
value 2. 

An alternative line of interpretation of the ob- 
served energy scaling laws could involve multi- 
scale Alfvenic wave packets excited in the chro- 
mosphere and corona by the photospher i c foot 
point motions fe.g.. | Ofman et al.l Il998t lUlrich 



1996 [ iKudoh k Shibata!ll999t ISuzuki k Inutsuka 
2006 :lDe Pontieu et alj2007llFuiimura k Tsuneta 



20091: lOfmanl 120101) . Perturbations in the photo- 
sphere are able to generate a broad spectrum of 
transverse waves with the statistical properties 
of the generating convective events. Once these 
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waves arrive at the corona they should result in 
dissipation events with similar statistics. 

In the Wentzel-Kramers-Brillouin approxima- 
tion, the energy flux of an Alfven wave packet gen- 
erated by a single shuffling event associated with 
the displacement velocity Vd is given by 

SW=pv*V A (34) 

in which p is the ambient mass density and Va 
is the local Alfven speed. Assuming that the dis- 
placement velocity Vd ~ (vh) is roughly constant 
(see the Brownian walk argument above) and the 
plasma is homogeneous, integration of (|34|) over 
the emitting area and time yields the Alfven en- 
ergy 

W A oc (v h f (B n ) [ dt [ dxdy (35) 

Since (B n ) ~ E/V and the domain of the integra- 
tion is defined the swept area S, we arrive at the 
scaling law 

W A oc l De - Dv+Ds = L Dw (36) 

which is identical to (|30|) . Thus, the wave mecha- 
nism leads to the same geometric scaling of the 
energy input from the photosphere as the one 
based on the nonpotential MHD Poynting flux. 
Both mechanisms predict D\y rj De and are in 
an agreement with the statistics of coronal en- 
ergy dissipation presented in section [2] The wave 
mechanism is it more likely to play a role in the 
generation of short-living coronal events while the 
long-lasting events could be driven by slow con- 
vective shuffling. Understanding the detailed in- 
terplay between these effects in the context of the 
multiscale photosphere - corona coupling is an im- 
portant task for future research. 

4. Conclusions 

The presented statistical results reveal a strik- 
ing similarity in the behavior of coronal emis- 
sion events (as measured by transient localized 
enhancements of STEREO EUVI emission flux) 
and photospheric magnetic events (as represented 
by the enhancements of SOHO MDI line of sight 
magnetic flux), when the populations in the two 
solar regions are matched by appropriate thresh- 
olds. We found that the event measures integrated 



over space and time (the swept area S, the spa- 
tiotemporal volume V, and the integrated flux E) 
arc described by essentially the same probabil- 
ity distribution functions, yielding the same val- 
ues of the power- law exponents rg, ry, and te 
for the photosphere and corona. The event mea- 
sures that do not involve spatiotemporal integra- 
tion show greater statistical discrepancies result- 
ing in differing log-log slopes for the two solar re- 
gions. Likewise, the geometric exponents describ- 
ing power-law relations between the measured pa- 
rameters and the linear size L of the event also 
have distinct values. The analysis of scaling re- 
lations between the obtained scaling indices has 
shown that the probabilistic behavior of the MDI 
and EUVI events is consistent with their charac- 
teristic spatiotemporal geometries and indicates 
a strong physical interaction between the photo- 
sphere and corona. The derived scaling laws are 
in an agreement with the classical model in which 
the photospheric flux tubes undergo quasi one- 
dimensional horizontal convcctivc motions driven 
by uncorrelated, random walk - like displacements 
which energize the corona in a multiscale fash- 
ion resulting in the power-law distribution of solar 
flares. 

Until now, the multiscale behavior of coronal 
dissipation has been considered in frames of two 
alternative statistical-physical frameworks. The 
first framework, self-organized criticality, seeks to 
explain power-law distributions of flare parame- 
ters via cooperative interactions of a large num- 
ber of nonlinearly coupled degrees of freedom 
representing unstable coronal loops and/or loop 
strands. If the interaction time scale is signif- 
icantly shorter than the time scale characteriz- 
ing photospheric motions, the interaction leads 
to scale-free avalanches of energy dissipation in- 
dicating that the system resides in the vicinity of 
a global critical state. The SOC behavior is ro- 
bust with respect to the changes in microscopic 
physics: the same scaling laws are expected for dif- 
ferent instability mechanisms as long as the critical 
dynamics remains in the same universality class 
(jSethna et al.ll200ll ). This property of SOC model 
can be invoked to explain statistical universality 
of different classes of flares. 

The second well-developed concept is the high- 
Reynolds number fluid turbulence. It has been 
proposed that multiscale energy release in the 
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corona could result from a turbulent intcrmit- 
tency (strongly inhomogeneous energy dissipa- 
tion) associated with energy cascade in t he corona 
(|Boffetta et alJll999t iNigro et all 12004) . In this 
framework, energy injection at the largest scale 
is due to photospheric motions whereas nonlinear 
effects in the corona are responsible for transfer- 
ring the energy to ever smaller scales, similarly to 
the generation of multiscale eddies in turbulent 
fluids. Up to the accuracy of the available instru- 
mentation, both SOC and turbulence are in an 
agreement with the observed scaling laws hinting 
at a possibility of a m ore universal phenomenology 
(|Uritskv et al.ll2007l ). 

An important underlying assumption shared by 
the two concepts is that the photosphere plays a 
role of a "structureless" driver, with the observed 
hierarchy of scales generated self-consistently at 
the coronal altitudes. Both groups of models are 
largely insensitive to the driving conditions as long 
as the scale separation requirements mentioned 
above are met. 

The results of this study lead to a drastically 
different interpretation which seems to be much 
simpler. We presented a case suggesting that it 
is the multiscale photospheric motion, and not a 
complicated intrinsic process in the corona such 
as SOC or turbulence, that shapes the statistical 
distribution of solar flares, at least for a quiet Sun. 
In other words, the coronal flaring activity has a 
multiscale spatiotcmporal pattern because the free 
energy input from the photosphere is fragmented 
into multiscale packets with a matching power-law 
probability distribution of sizes. The Occam's ra- 
zor principle requires that among competing hy- 
potheses, it is justified to pick the one which makes 
the fewest assumptions. Following this principle, 
the direct interaction between the two solar re- 
gions should be preferred. 

The proposed mechanism of multiscale proba- 
bilistic coupling is in a good quantitative agree- 
ment with scaling properties of the estimated free 
magnetic energy supplied to the corona by the 
photospheric shuffling and/or through an excita- 
tion of multiscale Alfven wave packets as shown 
in section 13.41 Each successful shuffling event 
generating a positive line of sight Poynting flux 
sends a packet of free energy into the corona via 
a multitude of DC and/q r AC co upling mech- 
anisms (see e.g. Davilal ( 1987 ): Ofman et al 



(Il994h: lOfman fc Davilal (fl995h ; iKlimchukl (|2006h : 
Mcintosh et all (|201lh ). According to our analy- 



sis, the amount of the injected free energy follows 
a power-law distribution described by the same 
log-log slope as the distribution of the integrated 
emission flux in the corona. It is worth mention- 
ing that the wave coupling scenario can reach far 
beyond the corona and directly affect the a ccelera- 
tion a nd heating of the solar wind (see e.g. lOfman 
l201dh . In this context, it would be important to 
investigate in future studies how the multiscale 
structure and dynamics of the solar photosphere 
are imprinted in the solar wind turbulence. 

It should be also noted that our conclusions ap- 
ply to relatively strong coronal dissipation events 
due to the thresholding which eliminated most of 
the small-scale flares. Generally, any threshold- 
based technique introduces an ambiguity with re- 
gards to the sizes of the detected features: the 
lower the threshold, the larger is the apparent size. 
The threshold dependence is especially strong 
when the probability density function of the mea- 
sured parameter exhibits a slow non-exponential 
decay as is the case with the distributions of MDI 
and EUVI pixel counts (see Fig. Since the 
events were selected by applying rather high in- 
tensity thresholds tuned to intermittent spikes in 
the studied images, the actual physical dimen- 
sions of the events in our analysis could signif- 
icantly exceed the estimated values. Therefore, 
our statistics are likely to describe relatively large 
events produced by strong magnetic discontinu- 
ities and involving considerable plasma volumes, 
possibly on multiple interwoven coronal loops. For 
smaller events occurring inside individual loops, 
the photospheric length scales are obviously not be 
directly "transferred" because the magnetic field 
strength decreases by orders of magnitude from 
the photosphere to the corona, and the transverse 
scale expands accordingly. The multiscale dis- 
sipation pattern appearing at these scales must 
build up in the corona independently of the pho- 
tospheric structure, e.g. following SOC and/or 
turbulence scenarios. 

In summary, the conducted analysis strongly 
suggests that energy dissipation scales of coro- 
nal heating events can be directly controlled by 
turbulent photospheric convection. Our findings 
speak in fa vor of the coupling scenario proposed 
by Parker ( Parker 1983L 1988 ) in which random 
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photospheric shufflings generate marginally stable 
magnetic discontinuities at the coronal level. They 
are also consiste nt with the the Alfven wave heat- 
ing mechanism dDavila 1987 ; Ofman et al J 1 19981 : 
De Pontieu et al.l 12007 ) which can play a major 
part at shorter time scales. Subsequently coro- 
nal instabilities relax a subset of unstable mag- 
netic structures or Alfvenic wave packets repeat- 
ing their probabilistic patterns. The value of the 
energy distribution exponent (je ~ 1-5) result- 
ing from this coupling is too low for the smallest 
dissipation events to play a prevailing role in the 
coronal heating process. The large-scale processes 
such as those driven by the mesogranulation flow 
are likely to be releasing a dominant portion of 
free magnetic energy dissipated in the corona, at 
least for the range of coronal emission intensities 
considered in our study. The question remains: is 
this energy sufficient to provide the bulk coronal 
heating necessary for the observed radiative loss 
rates? 
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Table 1: Percentile thresholds used to detect intermittent events in SOHO MDI and STEREO EUVI image 
sets (n - number of detected events). 



p 


MDI threshold 


EUVI threshold 


95.0 % 


13.5 (n 


= 38242) 


202 (n = 


4124) 


97.0 % 


17.3 (n 


= 19114) 


217 (n = 


3005) 


99.0 % 


32.3 (n 


= 5912) 


258 (n = 


1269) 


99.5 % 


42.0 (n 


= 3410) 


292 (n = 


= 686) 



Table 2: Scaling exponents of the detected photospheric and coronal events, as a function of the percentile 
threshold. 





MDI 


EUVI 


X 


n 




A 


n 


T~x 


A 

^x 


p = 


95% 














L 




1.00 ± 0.00 


3.93 ± 0.01 


2.94 


1.00 ± 0.00 


2.68 ± 0.22 


1.66 


T 




2.02 ± 0.17 


2.31 ± 0.07 


2.64 


1.33 ± 0.09 


2.12 ± 0.11 


1.50 


A 




1.42 ± 0.03 


2.99 ± 0.06 


2.82 


1.59 ± 0.03 


2.03 ± 0.10 


1.64 


S 




2.74 ± 0.16 


2.06 ± 0.03 


2.92 


1.99 ± 0.04 


1.86 ± 0.05 


1.71 


V 




3.18 ± 0.19 


1.92 ± 0.02 


2.93 


2.58 ± 0.09 


1.64 ± 0.04 


1.66 


E 




3.65 ± 0.20 


1.80 ± 0.02 


2.91 


3.04 ± 0.05 


1.50 ± 0.04 


1.53 


P = 


97% 














L 




1.00 ± 0.00 


3.67 ± 0.19 


2.67 


1.00 ± 0.00 


2.45 ± 0.15 


1.44 


T 




2.15 ± 0.09 


2.00 ± 0.03 


2.14* 


1.38 ± 0.06 


1.92 ± 0.09 


1.27* 


A 




1.40 ± 0.02 


2.80 ± 0.02 


2.52 


1.60 ± 0.03 


1.82 ± 0.02 


1.31 


S 




2.82 ± 0.13 


1.91 ± 0.03 


2.56 


2.00 ± 0.03 


1.80 ± 0.01 


1.59* 


V 




3.32 ± 0.14 


1.75 ± 0.03 


2.49 


2.61 ± 0.04 


1.60 ± 0.04 


1.57 


E 




3.77 ± 0.12 


1.65 ± 0.03 


2.46 


3.07 ± 0.03 


1.43 ± 0.04 


1.31 


P = 


99% 














L 




1.00 ± 0.00 


3.08 ± 0.03 


2.12 


1.00 ± 0.00 


2.11 ± 0.17 


1.08 


T 




2.21 ± 0.06 


1.84 ± 0.04 


1.84* 


1.34 ± 0.14 


2.04 ± 0.11 


1.39* 


A 




1.34 ± 0.04 


2.53 ± 0.03 


2.06 


1.60 ± 0.09 


1.77 ± 0.13 


1.23* 


S 




2.75 ± 0.09 


1.79 ± 0.04 


2.18 


1.99 ± 0.01 


1.75 ± 0.07 


1.50* 


V 




3.39 ± 0.10 


1.63 ± 0.04 


2.15 


2.63 ± 0.03 


1.56 ± 0.06 


1.47* 


E 




3.78 ± 0.10 


1.52 ± 0.04 


1.96 


3.01 ± 0.02 


1.47 ± 0.06 


1.42* 


P = 


99.5% 














L 




1.00 ± 0.00 


3.22 ± 0.35 


2.30 


1.00 ± 0.00 


2.53 ± 0.01 


1.54 


T 




2.14 ± 0.08 


1.74 ± 0.03 


1.58* 


1.28 ± 0.04 


2.32 ± 0.14 


1.69 


A 




1.43 ± 0.01 


2.50 ± 0.10 


2.14 


1.62 ± 0.11 


1.78 ± 0.15 


1.26* 


S 




2.81 ± 0.19 


1.68 ± 0.03 


1.90* 


1.95 ± 0.03 


1.86 ± 0.10 


1.68 


V 




3.46 ± 0.21 


1.49 ± 0.07 


1.68* 


2.60 ± 0.00 


1.58 ± 0.07 


1.51 


E 




3.82 ± 0.22 


1.44 ± 0.04 


1.68* 


2.97 ± 0.06 


1.58 ± 0.10 


1.71 
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Tabic 3: Average values of scaling exponents obtained by combining p = 99% and 99.5% data for the 
photospheric events, p = 95% and 97% data for the coronal events. The distribution exponents of time- 
integrated parameters are highlighted with bold font. 



X 


D, 




D, 






T x 










t x , predicted 


L 


1.00 


± 


0.00 


1.00 


± 


0.00 


3.15 


± 


0.18 


2.57 


± 


0.13 


2.55 


T 


2.18 


± 


0.05 


1.36 


± 


0.05 


1.79 


± 


0.03 


2.02 


± 


0.07 


1.91 


A 


1.39 


± 


0.02 


1.60 


± 


0.02 


2.52 


± 


0.05 


1.93 


± 


0.06 


1.95 


S 


2.78 


± 


0.11 


2.00 


± 


0.03 


1.74 


± 


0.03 


1.83 


± 


0.03 


1.74 


V 


3.43 


± 


0.12 


2.60 


± 


0.05 


1.56 


± 


0.04 


1.62 


± 


0.03 


1.56 


E 


3.80 


± 


0.12 


3.06 


± 


0.03 


1.48 


± 


0.03 


1.47 


± 


0.03 


1.42 
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Fig. 2. — The first (top) and the last (bottom) images in the studied conjugate sets of SOHO MDI and 
STEREO EUVI observations (left and right columns, correspondingly). 
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Fig. 3. — Probability distributions of MDI magnetic fluxes (top) and EUVI emission fluxes (bottom) describ- 
ing the entire observational interval. Solid, dotted, dashed and dot-dashed vertical lines show respectively 
the percentile levels p =95.0, 97.0, 99.0, and 99.5 % used to detect intermittent events in each data set.). 
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Any L 3 Mm< L < 12 Mm 




• MDI EUVI • MDI • EUVI 

Fig. 4. — Spatiotcmporal plots showing about 1/4 of the photosphcric (MDI) and coronal (EUVI) events 
detected at respectively p = 99% and 95% percentile levels, (a) No restriction on the liner size L of the 
events is applied; (b) linear size is limitted by the interval of scales L £ (3, 12) Mm used for evaluating 
scaling exponents. 
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Fig. 5. — Scatterplots of photosphcric magnetic enhancement events (left) and coronal emission events (right) 
detected respectively from SOHO MDI and STEREO EUVI images at the 95% percentile threshold. 
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Fig. 6. — Probability distributions of the photospheric magnetic events (left) and the coronal emission 
events (right) detected respectively from SOHO MDI and STEREO EUVI images at four different percentile 
thresholds shown in Fig. [3] 



28 




E (arb. units) E (arb. units) 

Fig. 7. — Comparison of probability distributions of photospheric and coronal events obtained for two 
combinations of thresholds yielding comparable numbers of detected events in each data sets. The slope 
values shown on each panel are the average log-log slopes for each pair of MDI and EUVI distribution 
evaluated over a combined range of scales defined by a union of the MDI and EUVI scaling ranges shown in 
Fig. [H^a) and (b) with vertical lines. 
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Fig. 8. — Temporal evolution of power-law probability distribution exponents of MDI and EUVI events 
computed within a 90-minute moving window. Vertical bars show the standard errors of the plotted power- 
law indices. 
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Fig. 9. — Schematic illustration of the characteristic spatiotemporal geometry of magnetic and emission 
events consistent with the conducted scaling analysis. Labels (1) and (2) refer respectively to the initial and 
final locations of the event; the shaded region is the projection of the event onto the image plane defining 
the swept area S. 
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